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Abstract 



Contributions to HFS and to the Lamb shift intervals of order a 2 (Za) 5 
induced by the graph with the two-loop overlapping electron self-energy dia- 
gram inserted in the electron line are considered. Explicit expression for the 
overlapping two-loop self-energy diagram in the Fried- Yennie gauge is ob- 
tained. Contributions both to HFS and Lamb shift induced by the diagram 
containing such subgraph are calculated. 



1 Introduction 



Calculation of contributions of order a 2 (Za) 5 m to hyperfine splitting (HFS) 
and Lamb shift induced by the two-loop radiative photon insertions in the 
electron line was initiated in the previous paper We have calculated 
there the contributions induced by all graphs containing one-loop electron 
self-energy diagram as a subgraph, by the graph containing two one-loop 
vertices, and also the contribution induced by the rainbow two-loop electron 
self-energy insertion in the electron lineQ. Calculation of the contribution 
induced by the overlapping two- loop electron self-energy insertion (see Fig.l) 
is presented below. Consideration of this contribution is impeded by all 
the usual difficulties which are connected with the presence of overlapping 
infinities in this case and the necessity to perform all necessary subtractions. 
In Section |6] we obtain explicit expressions for the renormalized overlapping 
diagram contribution to the self-energy operator in the Fried- Yennie (FY) 
gauge in the form of the Feynman parameter integral. These expressions 
are used later for calculation of the contributions to HFS and Lamb shift. It 
should be mentioned that taking into account also expressions for the rainbow 
and one-particle reducible [||] diagrams we have explicit general expression 
for the two-loop mass-operator in the FY gauge. We are aware of only one 
other calculation of the two-loop mass operator existing in the literature 
[[| which was performed in the Feynman gauge. We hope that the explicit 
expression in the FY gauge presented below may find applications also in 
other problems, e.g. in calculation of the positronium decay rate [|J. 



1 We use a chance to correct an error made in M in the calculation of the rainbow 
diagram contribution to the Lamb shift. Due to a normalization error the result cited in 
|l[ is two times larger than the correct one 



We are deeply grateful to Dr. Pachucki who attracted our attention to this error. 
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2 General Expression for the Overlapping Self- 
Energy Operator 

Consider contribution to the two-loop self-energy operator induced by the 
diagram with overlapping photons. First of all one has to perform subtrac- 
tions of potential infinities in internal vertices. Explicit expression for such 
internally subtracted contribution to the mass-operator in the FY gauge has 
the form 

{ p + l + m p + l + q + m p + q + m 
7a| D(p + I) 7/3 D{p + I + q) 7A Dip + q) 

p + I + m po + l + m p + q + m 



D(p + l) lf3 D(p + l) 1X ~D(p + q) 
p + I + m p + q + m p + q + m 



£>(p + ,p ^( Po + g) ^(Po + g) 



where 



= p 2 — m 2 



and po — m is the mass-shell momentum in subtracted terms. 

There is a certain subtlety connected with this mass-shell limit. Although 
renormalization procedure in the FY gauge is performed below without in- 
troduction of the infrared photon mass (see, for more details 0]) there are 
spurious infrared divergences on the intermediate stages of calculations on 
the mass-shell which may produce finite but discontinuous on the mass shell 
results if the mass-shell limit is taken in the naive way prior to calculation of 
the integrals. To avoid these problems one has to perform calculation even 
of the subtracted terms with slightly off mass-shell external momentum (in 
this case with off mass-shell momentum po) and only after calculation of all 
infrared unsafe integrals to take the mass-shell limit. In this way one pre- 
serves continuity of all physical results on the mass shell and obtains correct 
results. 
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We would like to obtain the expression for the overlapping diagram con- 
tribution to the self-energy operator in the form of the integral representation 
with minimal number of the Feynman parameters. To reduce the number of 
the parameters we transform the integrand in eq.fll]) with the help of trivial 
identity 



i - ,p + m ? .p + m 

p + l + q + m = (^—+l) + (^— + q) (2) 



and obtain 



1 p + q + m p + q + m 



where 



_ {p + 1 + m)n/p(p + 21 + m) _ (pg + j + m^ptpp + l + m) 
9 D{p + l)D{p + l + q) D*{p + l) ' U 

One may check that both terms in the braces in eq.(0) produce after in- 
tegration coinciding results, so we consider below only integral of the second 
term and simply double the result. Initial expression for further transforma- 
tions has the form 

W n / « n 2 f A 1 / fa , y^^n / ,p + q + m 



v 47r 7 iix 2 q 2 q 2 ' D(p + q) 

where 

/■ <i 4 Z 1 /°7 A 

^(p, ?) = J ~2p(9 aX + 2—) la n plx . (6) 

Integration over I in eq.(H) is convergent but separate terms in the inte- 
grand produce logarithmically divergent results. We slightly rearrange inte- 
grand in order to separate compensating divergences 

RpM= fU { l*D(p + l)D(p + l + q ) (7) 
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[(g aX + 2 l -^) la (p + 1 + m) 7 /j(p + 21 + m) lx - 6Z 2 7/3 ] 

[(g aX + 2—) la (p + I + m) lf3 (p + I + m) 7A - 3Z 2 7/3 ] 



Z 2 L> 2 (p + / 

1 1 
+67 ^ l D(p + Z)D(p + Z + g ) ~~ £> 2 (p + 

Representation in eq.(0) is very convenient for further transformations as 
each of the terms Rp (after substitution in eq.(^)) does not need more than 
four integration parameters to represent it in the form of the integral and 
each of these terms is a convergent integral in respect with integration over I 
which will be the first integral to perform. We will consider transformations 
of each of the terms in eq.(^) separately. 



3 Calculation of the Contribution to the Mass 
Operator Induced by the Term Iva 

Let us consider first contribution to the mass operator induced by term R$ 
in eq.(^). After ordering over powers of the integration momentum I we 
obtain 



. . . ( 8 ) 
+ [Qml p - 2ljpp - pjpl + 2l(p + 771)79 + 70 (p + m)l\ - [21^1 + Z 2 7/3 ]}. 

Next we combine denominators with the help of the identity 

(1 - x)l 2 + x[(l - z)D(p + /) + zD(p + l + q)] = (l + xQ) 2 - xA, (9) 
where 

Q = p + qz, 



4 



A = m — p (1 — x) — q z(l — xz) — 2pqz(l — x). 

After shift of the integration variable I — > I — xQ we perform momentum 
integration which leads to the expression 

Rf = 2 £ dx f a efe{2p / (2 -^- 3m + " (10) 



7/3 , ,QlpQ, , r x ^(i) , x 
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-(p» - m a )[(2 - + x(l - z)^] + + -Cf ]}. 



4 la) +< } +4 lc) . 



where 



Cf = 6mQ p - 2Q-f P p - fnpQ + 2Q(p + m) lp + 7/J (p + m)Q, (11) 

Cf = 2Q lp Q + Q 2 7/3 . 

Each power of momentum I in the numerator of the integrand in eq.(H) 
produces respective power of factor x in the numerator of the integrand in 
eq. (pTDj) . We will consider calculation of the terms in different brackets on the 
right hand side (i.e. terms Ra ) separately. 

3.1 Contribution to the Self-Energy Operator Induced 
by the Term 

It is easy to see that the expression in eq.(|T0|) leads to the contribution to 
the unrenormalized self-energy operator which remains finite even on the 
mass shell. However, we would like to obtain expression for the renormalized 
(doubly subtracted) self-energy operator and each subtraction inserts (after 
integration over q) additional power of x in the denominator, making inte- 
gration over x potentially unsafe. Only term in eq.(§) seems to produce 
infrared divergency on the mass shell after subtraction. To get rid of this 
apparent infrared divergency we first separate infrared unsafe terms in the 
expression for this term 

<» = 4 W p*/; ^ (2 -f- 3m +*(i-*)^i (12) 
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. /, ,QPQ-PQ 2 , 

+x(l-x) ^ }. 

Next we use trivial identity which is the result of integration by parts and 
the relation dA/dx = Q 2 

( dxxil — x)Q— — ( dx — - — (13) 
Jo A 2 Jo A v ' 



and obtain 



R p = A PP J dx J dz i 3 ^ + s(l - x) ^ }• (14) 

This new representation is completely infrared safe and admits subtrac- 
tion on the mass-shell. Respective contribution to the mass operator has the 
form 

Next we introduce Feynman parameters with the help of the identity 



(1 _ t )q 2 + t[u(- A ) + (1 - u)D(p + q)\ = (q + pqtf - M, (16) 
z{L — xz) 

where 



^ 99 / 9 9x/- \ m 2 — p 2 (l — x) m 2 (uj + p£) .„„. 

n = pWt + (m 2 - p 2 1 - u + / 1 , > u = — ^, 17 

z(l — xz) z(l — xz) 

x(l — z)u 

V= l ; , 

1 — xz 

1 1 

m — p 

P = 2 ' 

£ — (1 — rr)w + z(l — xz)(l — u — rft), 



6 



uj = ux + z(l — xz)r] 2 t. 
Denominators in eq. (|TB[) are combined as follows 



- . 1 — — - = f 1 du f 1 dt- gfeli^—— = G(k, l,n), (18) 

A k D(p + q)q 2n Jo Jo [(q + pqt) 2 - tQ} k+1+n V ' ' ; ' V ; 

where 



g(i, l, o) = - 1 , g(h l, l) = - n 2t y g(h l, 2) = - 6 * (1 % 

z(l — xz) z{l — xz) z(l — xz) 

(19) 

£2,1,1= — -, 2,1,2 =-27^ w- 

z z [l — xz) z Z Z {1 — XZ) Z 

We displayed in eq.(O) and eq.flTS]) slightly more general formulae than 
those which are necessary to perform momentum integration in eq. (|T5|) . 
These auxiliary equations will be extensively used below for momentum in- 
tegrations of other entries in eq.(||p. 

Contribution to the mass operator in eq.(|i~5D has a rather simple form in 
new notation 



E (ia) (p) = ( « f ['dx [ l dz[ ^| J2 N a (k,l,n)G(k,l,n), (20) 
4tt Jo Jo J m 2 k 

where 

A^(l,l,l) = I2\p(l-x)-m)[2p(p + m) + {2qp + pq)), (21) 

N a (l, 1, 2) = 12[p(l — x) — m](p + m)qpq, 

N a {2, 1, 1) = 4x(l - x)(Qj3Q - j3Q 2 )[2j3(j3 + m) + (2qp + pq)}, 

N a {2, 1, 2) = 4x(l - x){QpQ - pQ 2 ){p + m)qpq. 

We want to obtain after momentum integration as small powers of factor 
Q in the denominator as possible since in this case next integrations are 
more convergent and subtraction procedure is more accessible. Due to this 
reason we separated in the numerators characteristic structure QpQ — pQ 2 
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which does not contain cubic in external momentum p terms after shift of 
integration momentum and leads thus to smaller power of Q after integration. 

Next we shift integration momentum q — > q~pr]t and obtain after straight- 
forward integration in eq. (|T5|) 



(p) = (— ) 2 f 1 dx C dz f 1 du f 1 dt . 12 r {[p(l-x)-m] (22) 
Mvr 7o 7o Jo Jo z(1 - xz) 

{^[2p(p + m) - 3pV] - p(p + ™) + *(1 - 

+ ^— ^ {-^[p(l + 2«yt) - 2(p + m)z)- - t] (p + 3m) 

£> 4 fl — 

+ - V ^ 2 ; [2(p - m)^ - (3p - m) 2?? ¥]}}. 

One may get rid of denominator f2 2 with the help of integration by parts 
over t, taking into account that dQ/dt = p 2 if (compare eq. fll~3|) ) . Integration 
by parts leads to substitutions in eq . (Effi) 



t(l - t)p 2 r] 2 1 - 2t 
t 2 (l - t)p 2 r] 2 t(2 - 3t) 



(23) 



tt 2 tt ' 

Note that although the right hand side in eq.(^) is finite even on the 
mass shell, separate terms in the integrand may lead to infrared divergent 
contributions after subtraction. Integration by parts gave us the chance to 
get rid of these would be IR divergence which could become dangerous for 
the accessibility of the subtraction procedure. 

Now expression in eq.(E^) acquires the form 



ry rl rl rl rl 19 

^ la) (p) = (-) 2 dx dz du dt {3pt\p(l-x)-m\ (24) 

4tt Jo Jo Jo Jo m l [uj + p£) 

x(l — x)u 1 — 2t 
\p(l-T])+m]-\ — —p 2 \2{p-m) pt(l+2zr)t)+zt[m-5p(l-2t)}}}. 

1 — XZ 7] 

This form is very suitable for subtraction which we postpone until all 
other contributions to the mass operator would be obtained. 
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3.2 Calculation of the Contribution to the Mass Op- 
erator Induced by the Term Wp 

Contribution to the self-energy operator corresponding to has the form 
E( 16 >(p) = -2(f )V-™ 2 ) f dx f dz I '^f{(2 - xj% (25) 

47T JO JO J lit A 



A 2 q 2 D(p + q) q 2 

Again as in the previous section we try to transform numerator of the 
integrand to the form containing minimal power of external momentum p. 
This goal is easily reached with the help of trivial identities 



QlpQiP + q + m)^ = 2Q 2 (2p + 2q - m) + 2(1 - z) [QqQ - qQ 2 }, (26) 

QqQ = qQ 2 + [QqQ - qQ 2 ]. 

Note that the term QqQ — qQ 2 as well as the similar term QpQ — pQ 2 in 
the previous section is free of cubic in the external momentum terms after 
shift of integration variable. We get rid of the terms containing Q 2 in the 
numerator with the help of integration by parts displayed in eq. (|T3l) . We 
then obtain 



^)(p) = -2{^-) 2 {p 2 -m 2 ) f dx f dz [ ^ £ N b (k,l,n)G(k,l,n), 

47T JO JO J V7T k>n=1)2 

(27) 

where 

N b (l, 1, 1) = 6[g(l - 2x) -px + m(2 - a;)], (28) 

N b (l,l,2) = 6(l-x)gpq, 
N b {2, 1, 1) = 2x{l -x){2- z){QqQ - qQ 2 ), 

N b {2, 1, 2) = 2x(l - x){QqQ - qQ 2 ) (p + m)q. 

Performing next straightforward momentum integration one obtains 
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^m-^jM^jW* 2 ^ (29) 

{-\pr)t(l - 2x) + px - m(2 - a;)] + p(l - + *(1 - *)^|-] 

, ~ x)u ( p(2 - z)t _ p 2 (l-t) { p-m _ 

+ i-xz 1 n n 2 [ z mn ' 

Next we integrate by parts to get rid of denominators Q 2 (compare pre- 
vious section) and obtain 

S (16) (P) = (^) 2 tdx Cdz fdu f 1 dt{ 12 ^f ~ m 2 {P[x + (1 - x)(2 - 3t) 
47r Jo jo jo jo m z {uj + p4) 

(30) 

+?7t(l - 2x)] - m(2 - x) + ^0— — [p(2 - z)t + 2p l - 2t ]} 

1 — xz 1] 

_ ,2 ^P 2 (P + m)x(l - x)u{l - t) 

One term with denominator Q 2 ~ (a; + p£) 2 survived integration by parts 
but it already contains explicit factor {p — m) 2 and is completely safe under 
subtraction. 



3.3 Contribution to the Self-Energy Operator Induced 
by the Term 

Unlike previous terms the term is ultravioletly divergent and leads to 
a logarithmically divergent contribution to the self-energy operator 



(31) 

Since all divergent contributions contain extra powers of integration mo- 
mentum q in the numerator we transform numerators prior to integration 
separating powers of q 2 explicitly. Then we cancel this factor q 2 over similar 
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denominator factor reducing thus the number of Feynman parameters in the 
divergent terms. In the end of calculations we will insert additional integra- 
tion parameter in these terms in order to obtain integral representation for 
the divergent term which does not contain logs. As a result we will obtain 
integral representation for the contribution of the divergent terms to the self- 
energy operator which contains only four integration variables and contains 
denominators which are not too singular on the mass-shell . 

Extraction of the explicit factors q 2 is performed with the help of the 
identities 



C { p (p + m) (a p + = Cf (p + m)^ + 3p(-3p 2 + Amp + 2m 2 ) (32) 



+6z[m(qp — pq) + 3m(p + m)q — 2p q] 

+ (3p 2 + Imp + Am 2 ) 2 ^ + J q2 + 12*^^, 

q 1 q l 

\ 2qq^ 

Cp^q^P H — ) = q 2 3z(8m — hp) + 2(4m — 3p)(pq + 2qp) — z(2qpq + pq 2 ), 



Cf{p + m)(-f + ^f-) = Cf{p + m)^ + q 2 3z 2 (2m - p) + 3pp 2 



+2z(Ap 2 q+5mpq+mqp)+3z' 2 (2qpq+pq 2 )+p(p+2m) ' \-Az — , 

q 2 q l 

Cfqin P + ^f-) = q 2 (12z 2 q + 18zp) + 2z(2qpq + q 2 p) + 4(2pqp + p 2 q). 
Then we obtain contribution to the self-energy operator in the form 

^\p) = 2{^-) 2 ('dx fdz f^j2N c (l,l,n)G(l,l,n), (33) 

47T Jo Jo J t7T Z ~n 



a \2 f 1 j f 1 j f d 1 

i, 

where 



N c (l, 1,0) = 3xz[4qxz + p{-5 + 6x - xz) + 2m(4 + xz)}, (34) 



N c (l, 1, 1) = [xC { } ] + x 2 Cf ] ]{p + m)^ + 3px(-3p 2 + Amp + 2m 2 ) 
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+2x(4m — 3p){pq + 2qp) + 3pp 2 x 2 + Ax 2 {2pqp + p 2 q) + 6xz[m(qp — pq) 

+3m(p+m)q—2p 2 q]+2x 2 'zi^Ap 2 q+5mpq+mqp)+xz(—l+2x+3xz)(2qpq+pq 2 ), 

N c (l, 1,2) = [p 2 x(3+x)+2mpx(l+x)+4m 2 x](2qpq+pq 2 )+4xz(3+x)(pq)qpq. 

Performing next momentum integration in eq. (|3^) (note that in the term 
with numerator N c (l, 1,0) integration momentum is shifted as q — > q — pt] 
unlike the standard shift q —>■ q — pvjt in all other terms) we obtain 

£( lc )(p) = -2(-^) 2 C dx C dz f 1 du C dt—^ (35) 

^' v 4vr 7 Jo Jo Jo Jo z(l-xz) 

J/v/Vi 1 nV! A ' i^^2 2 A4(l,l,l) (l-t)AT c (l,l,2) 1 
{jVc(l,l,0)(bg^— -l) + 6px 2; 77 + — }, 

where 

A/" c (l, 1, 0) = 3xz[p(-5 + Qx-xz - Axzr]) + 2m(4 + :&z)], (36) 
W c (l, 1, 1) = pp 2 [9x(-3 + x)+ 6x{3 - 2x)r)t + 10xz(3 - 2x)r)t 
+3xz(— 1 + 2a; + 5a;2;)?7 2 t 2 ] + mp 2 6x[6 — (4 + 7z + 2xz)i]t] + m 2 p6:r(3 — Bzrjt), 
A/" c (l, 1, 2) = pp 2 r] 2 t 2 [p 2 x(3 + x)(3- Azrjt) + Qmpx(l + x) + 12m 2 x] 
and 

n lt=1 = n(t = i). (37) 

Once again one may easily get rid of denominator Q 2 with the help of 
integration by parts over t. We then obtain 



TS Xc \ P ) = -(^) 2 2 dx ['dz I 1 du \\t{ M f^ X % \og^- - 1) (38) 
KFJ v 4tt ; Jo Jo Jo Jo l z(1 - xz) v tt lt=1 



Z(1-XZ) li| t=1 

6pxVr? W c '(l,l,2) -A/" c (l,l,l) 



z(l — xz) m?{uj + p£) 

where 



A^(l,l,2) = 3p(l-2t)a;[p 2 (3+a;)+2mp(l+a;)+4m 2 ]-4pp 2 (2-3t)x(3+x)^t. 

(39) 
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4 Calculation of the Contribution to the Mass 

(2) 

Operator Induced by the Term R\ 

Contribution to the mass operator connected with the term R% in eq. (0) is 
the simplest one to calculate. Let us start with convergent integration over I 

+ {Aml p + 8(p Z)7 j8 ] - [2/V + l 2 1/3 ]}. 

We combine denominators in the same way as in eq.([|) (z = now since 
there are only two different factors in the denominator) 

(1 - x)l 2 + xD(p + /) = (/ + xp ) 2 - a;Ao, (41) 
where A = A(z = 0). 



Integrand on the right hand side in eq.(40) is a sum of terms contain- 
ing increasing powers of integration momentum I (or powers of the Feyn- 
man parameter x after integration over momentum). Here we encounter the 
problem concerning the term with the lowest power of integration momen- 
tum. It is easy to see that if one puts momentum p to be exactly on the 
mass-shell then respective integral turns out to be the sum of two infrared 
divergent terms. These divergences, of course, cancel leading to a finite re- 
sult, however, one has to be extremely careful performing this calculation. 
Correct way to perform integration here is to preserve small nonzero virtu- 
ality po — ( m2 — Po) I m2 connected with momentum po in the denominator 
A = m 2 [x + po(l — x)] at intermediate stages of integration. We thus reg- 
ularize would be infrared divergences, perform then necessary integrations 
and go on the mass-shell only in the end of all calculations 



Rf = -2m 2 7 , [ dx{-A[^ - 2m 2 ^-^} + 12^ + 3^} (42) 

JO ^0 ^0 ^0 

= -37/3- 

This simple integration nicely illustrates the necessity to be extremely 
careful about going to the mass-shell limit in the presence of the would be 
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infrared divergences. Really, if we naively put po = (m, 0) prior to integration 
over x in eq. (f42|) we would obtain a finite result (potentially infrared divergent 
terms cancel in the integrand even in this case) but continuity in the external 
momentum would be lost and, hence, that result would be wrong. 

It is a simple exercise now to obtain respective contribution to the mass 
operator 

It is easy to see that integral in eq.(|43"D only by numerical factor differs 
from the expression for the unrenormalized one-loop self-energy operator in 
the FY gauge (see, e.g [|]). Hence, we may use results of one-loop calcula- 
tions and immediately put down contribution to the renormalized self-energy 
operator induced by the term under consideration 

^fr) = (44) 
and we remind that virtuality p was defined in eq . (|T7|) . 



5 Calculation of the Contribution to the Mass 

(3) 

Operator Induced by the Term Rg J 

(3) 

Consider contribution to the mass operator induced by the term Ra in 
eq.(0). Integration over I is again convergent and one easily obtains 

(3) f d^l 1 1 

Rp =&1P i ^ [ D(p + l)D(p + l + q) ~ D 2 (p + l)^ (45) 

f 1 , _ m 2 f 1 z(l — 2z) 

= 6"f0 / dzlog — = -67/3 / dz , 

Jo Ai Jo Ai 

where 

Ai = A(x = 1) = m 2 - q 2 z(l - z). (46) 
Respective contribution to the self-energy operator has now the form 
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= -6(^) 2 dzz{\ - 2z) [ ^§{3(2m - p)G 1 (l, 1, 0) 

AlX JO J ITT 

where subscript accompanying function G\(l, 1, 0) means that one has to 
substitute in the respective definition in eq.fllBD Ai from eq. (|46|) . 
One may easily perform momentum integration now and obtain 



S (3) (P) = 18(f) 2 fdz 1 —^- fdu f dt{(2m - p)(log^-^ (48) 
An Jo 1 — z Jo Jo siii 



A 2 pp 2 {l -uft 2 



A 2 + iln 0|a, = i 

where 



tt lx=1 = n(x = 1) = p 2 (l - uft + (m 2 - p 2 )(l - it) + (49) 

2 

Q u =n(x = l,t=l)= p 2 (l - m) 2 + (m 2 - p 2 )(l — it) + — r. 

2(1 -z) 

Note that one has to preserve function fi n even on the background of 
the infinite cutoff A 2 to facilitate convergence of integration over Feynman 
parameters in eq.([4S|)p|. 

6 Total Expression for the Overlapping Dia- 
gram Contribution to the Self-Energy Op- 
erator 

2 Compare detailed discussion on this point in pj. 
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6.1 General Representation of the Overlapping Dia- 
gram Contribution to the Self-Energy Operator 

Next task is to perform double subtraction in all contributions to the self- 

). Con- 



energy operator obtained above in eq. fl24D , eq . (p0|) , eq. (p8f) and eq.(| 
tribution in eq. (PHI) is already subtracted and does not need any further trans- 



formations. The term containing factor (p — m) 2 in eq.(|30"D is also already 
subtracted. All other contributions to the self-energy operator either contain 
in the integrand denominator uo + p£ or logarithm of the cutoff momentum. 
Consider first logarithmically divergent contributions. 
Divergent contribution in eq.(p8|) has the form 



4^(p) = -(f) 2 / dxf dzf du-^[(p-m)h 1 + mh 2 }(log- 1), 

47T JO J J 1 — XZ "|t=l 

(50) 



where 



hi = — 5 + 6x — xz — 4xzr], h 2 = 3 + 6x + xz — 4xzr]. (51) 

Momentum cutoff disappears after subtraction and emerging logarithm 
may be put down with the help of integral representation 

log^ ! — =f>€\t=i : — 7 • 52 

u\t=x Jo uj\ t =i + pt,\t=\v 

This parametric representation of the logarithm gives one a chance to 
gain an additional apparent factor p — m in the numerator which is very 
convenient for the subtraction procedure. 

Subtracted expression in eq.(pOf) acquires then the form 



«W = -(f ftefdzfiu^ (53) 

47T m Jo Jo Jo 1 — xz 

f 1 — c \p(h + 2h 2 ^) + m{hi + h 2 + 2h 2 ^)). 

Jo uj\ t =i + p$\t=iv u\t=i uj\ t=1 

Consider next divergent contribution in eq. (ffgp . It looks like exactly as 
the previous divergent contribution with h\ = 1 and h 2 = —1 
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Sffi(p) = -18(^) 2 /* dz f du[(p -m)-m\ (54) 

47T JO 1 — 2 Jo 

A 2 + fl n _ A 2 
l ° g fi n A 2 + fi n j ' 
One easily obtains doubly subtracted contribution to the mass operator 



^ra / x / ol x0 (p — m) 2 f 1 , 1 — 2z r l , f 1 . . 

S^ fl (p) = -18 -H „ J / tfe- / rf<n / —— 55 

- 2^) - 2m^ 



where 



£n =((x = l,t = l) = z{\ - z)u{\ - u), 

lo\\ = uj(x = l,t = 1) = u + z(l — z)(l — u) 2 . 
Now we have to consider subtraction of all other entries in eq. ( plf ) , eq. (|30[) . 



eq. (|38|) and eq.(|48l) which contain factor uj+pt; in the denominator. Collecting 



all such terms we obtain the expression 



q; rl rl rl rl ]_2 (p — TTl\ 2 
S |^+P?(p) = ("J-) / rf;E / ^ / rfu / rft 1 5 [p/lp + Wl/lm] 

47T Jo Jo Jo Jo oj + p£ 



+(P - m )f2 + mf 3 }, 



(56) 



where 



/{ p = 2 + a;(-7 + x + 3t + xt) + [-2 + a;(4-2x + 9z-2xz-6zt-2a;2t)]^ + 5, 
f[ m = 2+ x (-8+x+5t+Axt) + [-l+x(A-Ax+llz-6xz-12zt-Axzt)}7]t+2S, 
f 2 = x (-7 + x + 8t + 7xt) + [-1 + x(Q - Qx + 8z - Wxz - I8zt - 6xzt))rjt + 3S, 
/a = 3x(l + z)t + ac(2 - 2x - 3z - Axz - 6zt - 2xzt)rjt + S, 
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f" 


= t(l 


f" 

J lm 


= t(2 - 




= f(l- 




tt = 




S = 



1 -2t 

V 

1 - 2t 

6 , 

rj 

,l-2t 
> , 

Tj 



2 v 

After subtractions we obtain 



4n m z Jo Jo Jo Jo uj + pt, uj 



+2/ 3 (-) 2 ] + m[f lm + (f 2 + / 3 )- + 2/ 3 (-) 2 ]}. 



(58) 



Subtraction in the last term in eq. fl48|) may be easily performed finishing 
thus calculation of all contributions to the subtracted self-energy operator 
induced by the diagram with overlapping photons. 

Collecting all entries to the self-energy we obtain 

e*(p) = (^) 2 (p ~T )2 f dx t dz t du t dt t d v{^—Av\hv + h- 

4tt m 2 Jo Jo Jo Jo Jo uj + p£ to 

(59) 

+2/ 3 (i) 2 ] + ,n[f lm + ( h + /,)« + 2/ 3 (i) 2 | - ^ + <"Ml-y-«) } 

uj uj m z (uj + pt,) 



1 - XZUJ\ t=1 + p£\ t=1 V U\t=l ' W\t=l 

_ lg z(l - 2z)(l - U )H" ^ + 3 ^ + 2( ^=i )2] + 2m(1 + fci_ )2} 

0J\x=l + P£|x=l ^[x=1 ^|x=1 ^|x=1 

61 .[p(l-2^)-2m^]+9- ^ 1 



1-2 Wn + p^nf W11 W11 x + p(l - x) 
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= ^ — m 2 " — \pvp(p) + ma m {p)}. 

Note that despite the appearance each term on the right hand side is at 
most a four dimensional integral. 

6.2 Infrared Safe Representation of the Overlapping 
Diagram Contribution to the Self-Energy Opera- 
tor 

The representation for the overlapping diagram contribution in eq. (|59|) is 
quite suitable for calculation of the respective contribution to HFS but needs 
further transformations for the Lamb shift calculation. The problem emerges 
when one tries to calculate values of the subtraction constants cr p (0) and 
<7 m (0) which enter the expression for the Lamb shift (see below). Integrals 
containing terms f2^/^ and 2f^ 2 /uj 2 in eq.([5^|) diverge as logp when p goes 
to zero (we remind that according to eq.([H]) u = xu + zt(l — xz) + 0(xuzt)) 
and this divergence disappears only in the sum 

M(p) = f 1 dx f 1 dz f 1 du f 1 dt—^—\^- + 2(^) 2 ] (60) 
Jo Jo Jo Jo UJ + pt, OJ 0J 

which enters eq.(^). 

Throughout investigation of the integrand in eq.(|60|) with the help of the 
definitions in eq.(|l7|) and eq.(|57|) shows that only terms 

{h£)m = x (~ 7 + x ) z 0- - xz ) - (61) 
(M 2 ) IR = xt[(5 + x-3z- Axz)z 2 (l - xz) 2 + u 2 (5 - u)] 

lead to divergencies at vanishing virtuality. Each term in the subtracted 
numerator / 2 £ — (f2^)iR contains either factor xu or zt which suppress di- 
vergency at vanishing virtuality. In the same way each term in the other 
subtracted numerator / 3 £ 2 — (fsC, 2 )iR contains one of the factors (xu) 2 , (zt) 2 
or (xu)(zt) and respective integrals are also infrared safe at vanishing virtu- 
ality. 

We are going to transform integral representation in eq.(^) to a form 
containing integrands which lead only to the integrals finite at p = 0. Can- 
cellation of infrared divergent terms may be greatly facilitated by reducing 
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the multiplicity of integration. To reduce the number of integrations we note 
that infrared properties of the integrands does not change if one substitutes 
in the numerators in eq.floTp x — > duo/du and tu — > udu/dz. Using this 
fact we separate infrared divergent terms by substituting these derivatives in 
eq. (|6l|) and define new infrared safe functions 

M = M-V 2 u, (62) 

where 

d d 

V 2 = (-7 + x)z(l-xz)—-u—, (63) 

d d 
V 3 = t(5 + x — 3z — Axz)z 2 {\ — xz) 2 — — h xu 2 (5 — u) — . 

du dz 

Now M(p) in eq . (|60D acquires the form 

M(p)^M(p) + M 1 (p), (64) 

where M(p) differs from M(p) only due to substitutions / 2 £ — > f 2 ^ and 
f 3 ^ 2 — > /3^ 2 . Term M(p) is infrared safe and logarithmic infrared divergency 
is connected only with the term 

M 1 (p) = f 1 dx f 1 dz f 1 du f 1 dt—t—l— + 2 ^-} (65) 
Jo Jo Jo Jo w + p( W UT 



[ dx [ dz [ du [ dt [ dv[- — — + 2(1 — v 

Jo Jo Jo Jo Jo (u + pvt) 2 



T> 2 u 2T> 3 uj 



{u + pv£) 2 1 (u) + pv £) 3 ' 

We easily extract infrared safe part M^(p) from the term Mi(p) with the 
help of the trivial substitution 



UJ 



-pvZ+(co + pvt), (66) 



which is valid according to the definition in eq.(|T7|). 

First term on the right hand side in eq. (|66|) produces infrared safe con- 
tribution to M 1 (p) 
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ri ri ri ri ri v 2 i 

pM,( P ) = - P J o dxJ Q dzJ Q duJ Q dtJ Q dvvlj-^+Ail-v). 

(67) 

while the second term on the right hand side contains complete derivatives 
over one of the Feynman parameters which may be easily integrated. Hence, 
we obtain the following representation for the Mi(p) term 

M x {p) = pM^p) + M up {p) + M down (p), (68) 

where 

M u Jp) = f 1 dx f 1 dz f 1 du f 1 dt f 1 dv{ ^~ X ^ l ~ XZ ^ (69) 
Jo Jo Jo Jo Jo u\u=i + pvq\ u =i 

u t(5 + x — 3z — Axz)z 2 (l — xz) 2 

— 2(1 — v)[- 



u\ z =i + pv£~\ z =i ~ (W|u=l + pv£\u=i) 

xu 2 (5 — u) 

(w| Z =i + pvi~\z=i) 2 



rl rl rl rl rl 

Mdown(p) = / dx dz du dt dv{ 
Jo Jo Jo Jo Jo 



-7 + x 



t + pv(l-t) 



1 +2(i-^)[ t(5 r +x - 3 ^-^ ) + r x{h 7 u) 



x + pv(l - x) [t + pv(l - t)] 2 [x + pv(l - x)Y 

It is easy to check that all integrands in M up (p) are infrared safe even 
at vanishing virtuality, while each integrand in the expression for the term 
Mdown(p) produces logp in the small p limit. Happily due to simplicity of 
parametric integrals in the expression for M down (p) one may easily obtain its 
parametric representation in the form of one-dimensional integral 



M, 



down 



(p) = ^/ 1 ^{-i-^lnl (70) 
KH) 2 Jo 1 2 1 - pv pv 



+2(1 - v)[-^- In 1 + _^_(^_ In 1 - 1)]}. 
1 — pv pv 1 — pv 1 — pv pv 

It is easy to see that total expression for M down (p) at p — > is free of logs 
and one may easily obtain 
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M down (0) = -H (71) 

Thus, we have obtained representation for the contribution of the over- 
lapping diagram to the self-energy operator in the form 



Sh(p) = fc) 2 ^— rH^(P) + mZ m (p) (72) 



a ^ 2 {p — m} 2 
'47r m 

+12(p + m)[pM^p) + M up (p) + M down (p)}}, 

where cr p (p) and a m (p) differ from the respective expressions in eq 
only by substitutions / 2 £ — > and / 3 £ 2 — > / 3 £ 2 (compare eq.(pl)). 



Each term on the right hand side in eq. (|72|) is explicitly finite at vanishing 
virtuality. Hence, this representation for the self-energy operator is very 
convenient for subtraction of its value on the mass shell and for calculation 
of the respective contribution to the Lamb shift. 



7 Contributions to the Energy Splittings in- 
duced by the Overlapping Mass Operator 

7.1 Contribution to HFS Splitting 

Radiative correction to hyperfme splitting induced by the graph in Fig.l is 
given by the matrix element of this diagram calculated with on mass-shell ex- 
ternal electron lines and projected on respective spin states and multiplied by 
the square in the origin of the Schrodinger-Coulomb wave function (for more 
details see, e.g. |]] an d @)- It is not difficult to obtain explicit expression 
to the energy splitting having formula for the mass operator in eq . ([59|)F1 



AE HFS = ^^-E F (-^) I " d|k| v P (k 2 ) (73) 
Tin 6 2n z Jo 



a 2 (Za) 1 u 9tt 2 



Tin 3 2ti 2 4 

where 



integration momentum k is measured in the units of electron mass in this section. 
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H= f 1 dx f 1 dz f 1 du f 1 dt{-^=[f lp + f 2 - (74) 

Jo Jo Jo JO v UJc UJ 



+2/3( i)2 ] + 3x(l-zMl-i),l 1 



UJ y/Uji £ UJ' 

2x 



.^L(3h 1 + 2h 2 ^L) 

1 - XZ\ U\ t =l U\ t =l 



-6^J^(3-2^)}. 
1 — 2; y Wil ^11 



Numerically we obtain 

AEhfs = -1-984(1) a ^ Z ^ E F . (75) 
7.2 Contribution to the Lamb Shift 

Contribution to the Lamb shift induced by the diagram in Fig.l is given by 
the matrix element similar to the one in the case of HFS, the only difference 
being in the spin structure and the necessity to perform an additional sub- 
traction on the mass shell, already mentioned above (see, e.g. [T] and 0). 
This subtraction is greatly facilitated by the representation in eq. ([72]) . Using 
this representation we obtain 



A£l = )M-f r *UM k2 ) + ? m (k 2 )] - MO) + g? m (0)]} 

7rn d 47r m Jo vl 

(76) 

+24k 2 M f (k 2 ) + 24[M up (k 2 ) - M up (0)} + 2A[M down (k 2 ) - M down (0)}} 



a 2 (Za) 5 ,m T .T_ 
-m' 



1 r°° 

— / rf|k|[A^ 1 (k 2 )+A£ L2 (k 2 )+A^3(k 2 )+A^ 4 (k 2 )]. 

i7T z JO 



irn 3 m 2i 
According to definitions in eq.(|59|) and eq.(|72|) 
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rl rl rl rl rl 

er„(k 2 ) + cr m (k 2 ) = / dx dz du dt dv{ 
Jo Jo Jo Jo Jo 
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rf f f t, ,te*f*?v 24(l-k 2 )x(l-xMl-t) 

^(1 -2z)(l -u)H 2 r „£\ x =i ,,61=1x2, 
-18— t~97 — — I 3 + 7 "^^ + 4(-^^) 2 ] 

W|x=l + K qa:=l <^|x=l ^|x=l 

_i 8 i^ ^ a _ + 9 - i 

l-z Lo n + k^ n v y con a; + k 2 (l -x) r ' 
Subtraction may be easily performed with the help of identities 

1,1 i, e 1 



±1 = -lr 1 + 1 

1,2 L/, . 1 l,2t\2 , ,2J , . L/, . 1 l,2t\2 1 



k 2 > + k 2 £) 2 cu 2J cu L (cu + k 2 2 cu(cu + k 2 
and we obtain 

rl p 1 ^1 ^1 ^1 

A^Li(k 2 ) = / (irr / dz I du I dt I dv 
Jo Jo Jo Jo Jo 

24a; (1 - x)u(l - t) r _ £ l9/ £ x9n 

+ - L pJ7 ! [2hi + h 2 + 4:h 2 - ) 

1-xz u\ t =i + k 2 £\t=iv u\ t =i ou\t=i 

+ i8 ^ 1 - 2 ^ )( 1 1 2 : M)2t2 fei[3 + rfci + 4(fci) 2 ] 

W|x=l + k 2 £|x=l W|x=l ^|x=l ^|x=l 

+l 8 i^ ^ ^(1 - 4^) - 9 — }. 

1 — z uJu + k 2 ^nw Ol>u Uu x + k 2 (l — x) 

Integration over v is greatly facilitated by the obvious relations 
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which are valid for arbitrary function F. After integration over u and k 
we obtain 



where 



9vr 2 

C /|k|A^ 1 (k 2 ) = -^ r + 7 rL, (81) 



r i ,1 ,1 ,i g ^ t 

L = / dx / tfe / d«/ =-Lfi„ + /i m + /3- (82 

jo Jo jo jo yoJt oo uj 



M + ^M!)i + 6x(i-xMi-t) + 

u; cu 2 cufv^ u; ' 



2x ^Wfc + ^V" 



1 - xz W|t=i 5 W| t= i 

^(1 — 2z)(l — «) 2 t 2 £|a; = i r „£b=l ,4=1,,, 
yJW\x=l£\x=l ^l* =1 W l*=l 

+6 lz^(|i)i (1 _^) } . 
1 — z 5 

Next we perform integration over v and |k| in the second term on the 
right hand side in eq.(|76|) 

roo poo 

/ d\k\AS L2 (k 2 ) = 24 / d\k\Mi(k 2 ) (83) 
JO jo 

= -24 f 1 dx f 1 dz f 1 du P dt C dvv [°° d\k\\- — V f n -r= 
Jo Jo Jo Jo Jo Jo 1 IL (cj + kX) 2 

+4(1 -v)-. = -4tt dx dz du dt—^[V 2 £+ ~ - Sl 



(tu + k 2 v£) 3 Jo Jo Jo Jo 5 tu 
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Contribution of the third term on the right hand side in eq.(^) is equal 

to 

! o °°d|k|A£ L 3(k 2 ) = 24 J™ ^[M up (k 2 ) - M up (0)) (84) 

"1 /•! rl rl rl rx v£\ u=1 {7 - x)z{l - Xz) 



rl pi rl rl rl roo 

24 / dx dz du dt dv d\k\{ 
Jo Jo Jo Jo Jo Jo 



o Jo Jo Jo Jo Jo ^\ u =i ^|tt=i + k 2 t>£| u=1 

vt^i u +2(1 _, )[ <h=i r 1 



to\ z =i uj\ z= i + k 2 <| 2=1 " W| u=1 (cj| u= i + k 2 <| u= 

' -]t(5 + £-3z-4:rz)z 2 (l-:rz) 2 



w(w|«=i + k 2 <| 



u=lj 



cu| 2= i (cj[ z =i + k 2 <| z= i) 2 lu(uj\ z=1 + k 2 v£,\ z=1 ) 



, , f 1 , /-i , 2(l-^)\/^«=i r „ 6t(5 + x-3z- 4xz)z(l -xz), 
8vr / da; / dz dt 3— ¥ {-7+x+ — i '-} 



o Jo Jo , .5 5 W| u= i 

Uy |M=l 



+8tt / dx / du / dt— V — {-1 + i -} 

Jo Jo Jo , ,2 5 wu=i 

= irL up . 

Last term on the right hand side in eq.flTB]) gives 

d|k|A£ L4 (k 2 ) = 24 r ^[M down (k 2 ) - M down (0)} (85) 
Jo Jo k z 



180 / rfw ^L{ln(k 2 u) - 2(1 - v)[\^v) + 1 + 

Jo Jo 1 — k^u 1 — k^t> 

/■°° ,„ , r Ink 2 4 r Ink 2 



12tt 2 . 



Combining expressions in eq. fl8"T|) , eq. (|53D , eq.(^|) and eq.(^) we obtain 
final expression for the contribution to the Lamb shift in the form 
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. ,_, ,m r .^a 2 (Za) 5 1 r 397r 2 , r 

l = ™(^) + n(L +L ^ + Lup)h (86) 

or numerically 

AE L = 1.749(2) ^f C^fm. (87) 
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Figure Caption 

Fig.l. Diagram with two external photons and overlapping two- loop electron 
self-energy insertion in the electron line. 
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Fig.l 
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